Block spin magnetism and metal-insulator transition in a two-dimensional Hubbard model with 

perfect vacancy superstructure 
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We study the phase diagram of a square lattice Hubbard model with a perfect vacancy superstructure. The 
model can be also defined on a new bipartite lattice with each building blocks consisting of a minimal square. 
The non-interacting model is exactly solved and a mid-band gap opens at the Fermi energy in the weak inter- 
block hopping regime. Increasing the Coulomb interaction will develop the Neel antiferromagnetic order with 
varying block spin moments. The metal-insulator transition with Umi smaller than the one without vacancies oc- 
curs above the magnetic instability Um- The emergent intermediate magnetic metal phase develops substantially 
in the moderate inter-block hopping regime. Drastic increases in the ordered moment and the gap magnitude 
are observed on the verge of tight-binding band insulator with increasing U. The implications of these results 
for the recent discovered (A,Tl)j,Fe2-i:Se2 compounds are discussed. 

PACS numbers: 



The discovery of superconducting iron pnictide materials 
have attracted enormous enthusiasm in searching for copper- 
free high temperature superconductors |{ll. Recently, a family 



of ternary iron chalcogenides (A,Tl)j,Fe2_2:Se2 Il2|,|3|] (with A 
being the alkali atoms)has been shown to exhibit supercon- 
ductivity (SC) in proximity to a magnetic insulating phase. In 
the iron deficient compounds (x > 0) there are Fe-vacancies 
in the square lattice of the Fe-atoms. Interestingly, the Fe- 
vacancies could be ordered in certain periodic superstructures 
in (K,Tl)Fe2-iSe2 where the SC appears for a; 0.12 ^ 0.3 
HQ]. The vacancy superstructures were reported in early 
Mossbauer experiment on the non-superconducting material 
TlFe2_a;Se2 some years agoyL and recent transmission elec- 
tron microscopy experiment H provides clear evidence for 
the Fe-vacancy superstructures in the KFe2_i^Se2 samples. 

First principle calculations for (K,Tl)Fei 5Se2 show that 
the Fe-vacancy orthorhombic superstructure can be stabilized 
in the ground states and play a ciTicial role in the electronic 
structures [7, 8]. The observed activation gap, which is about 
^ 60 meV for x ^ 0.5, raises a concern as whether the insu- 
lating behavior is due to the antiferromagnetic (AFM) order- 
ing itself, or due to the Fe-3d electron correlations. The band 
structure calculations in [8] provide evidence that a moder- 
ate Coulomb interaction is required in order to account for 
the sizable gap and a Mott localization can be enhanced due 
to the ordered vacancies. The paramagnetic (PM) Mott tran- 
sition for X = 0.5 as well as possible s-wave SC in a doped 
Mott-insulator have been also proposed in a two-orbital model 
using spin-rotor mean-field (MF) theorvfl fioll . 

Here, we focus on another type of Fe-vacancy pattern, 
i.e., a tetragonal vacancy superstructure shown in Fig. 1(a). 
The advantages of this pattern is of two-folds. First, it 
is a perfect vacancy superstructure in the sense that it ex- 
hibits the maximal symmetry without breaking the in-plane 
four-fold rotational invariance and that all sites are equally 
three-coordinated (with vacancy density 20%)[11J. Sec- 
ond, this structure has been clearly observed in a number of 



(K,Tl)j^Fe2-:rSe2 samples (for x ~ 0.4)10 [H [3. Recent 
first principle calculations lll ill show that the ground state of 
(K,Tl)j^Fei 6Se2 is a large block spin checkerboard or Neel 
AFM which is quite different to all previously known mag- 
netic patterns in the iron-based materials. This novel AFM or- 
dering pattern has been observed in recent neutron diffraction 
experiments 1. 1 3-1 5.1 . It is due to this large spin magnetism 
that the parent material (K,Tl)o.8Fei 6Se2 could be an AFM 
insulator lflTlfiell . 

These developments motivate us to study a simplified sin- 
gle orbital Hubbard lattice with the perfect vacancy super- 
structure. This model not only describes the Mott physics 
in the presence of the ordered vacancies, but also captures 
the geometry frustration as well as the structure distortion. 
Therefore, it can serve as a basic model system for the metal- 
insulator(MI) transitions and the complicated magnetic prop- 
erties in the presence of the vacancy orderings. In this paper, 
we show that the interplay among the electron correlation, the 
geometry frustration, and the structure distortion in such a va- 
cancy superstructure will result in a rich phase diagram. 





FIG. 1: (color online), (a) Square lattice with perfect vacancy super- 
structure: The dashed lines connect four neighboring vacancies en- 
closing a minimal square. All sites are equivalent with three coordi- 
nates. The X, y axes are for conventional Fe-Fe square lattice, while 
the X, Y axes are for the bipartite lattice of blocks, (b) The n.n. 
(red) and n.n.n. (blue) intra-block (solid) and inter-block (dashed) 
hopping amplitudes. 
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The Hamiltonian we studied is the Hubbard model 



H 



J2h 

ija 



E 



fj.{nif + n4),(l) 



where Cia- is the annihilation operator for a Fe-3d electron at 
site i with spin a =t,i; n-icr = cJo-Cio- the electron number 
operator; tij the hopping parameters; U the on-site Coulomb 
interaction; and fi the chemical potential. As a peculiar fea- 
ture of the vacancy superstructure, the whole lattice can be 
perfectly covered by fundamental square blocks enclosed by 
four neighboring vacancies under the periodic boundary con- 
dition, see in Fig. 1(a), where the X, Y are the new axes for 
the block lattice. Then the four sites in each block can be 
labeled hy A = (5m - 1/2, 5n - 3/2)/^/5, B = (5m + 
3/2, 5n - l/2)/%/5, C = (5m - 3/2, 5n + l/2)/V5, and 
D = (5m + 1/2, 5n + 3/2)/-\/5, where, integers (m, n) la- 
bels the block position and the shortest Fe-Fe distance is set 
to be unit. Hence the model is invariant under translations 
m — ?> m+lorrt — > n + 1 along either directions. We 
then let tij ~ ^1,^2 ( or tij ~ ^1,^2 ) be the intra-block ( 
or inter-block ) nearest neighbor (n.n.) and the next nearest 
neighbor (n.n.n.) hopping amplitudes, illustrated in Fig. 1(b). 
Notice that the n.n.n. hopping amplitudes are mainly due to 
the 3d-4p hybridization bridged by the Se-atoms above or be- 
low each blocks(not shown in Fig. 1(a)), and in general, the 
vacancies will lead to structure distortions |^,[lTl 16 1, break- 
ing the degeneracy between the intra- and inter-block hopping 
amplitudes. While these hopping parameters can be deter- 
mined by fitting the LDA band structures, our task here is to 
study the phase diagram for generic model parameters. 

To begin with, we first exactly solve the non-interacting 
tight-binding model, i/o = Eka. cL(^)^aa' (i?)cA.(i?), 
where, c\ry{K) is the Fourie transformation of cao-(p) for 
an electron at the site A in the p-th block, K is in the first 
Brillouin zone (BZ) of the periodic blocks. A, A' = A,B,C 
and D, and Hxx' = [0,ti + he"'^- M + iiC^'^y ,t2 + 



iie-'^y]; [h + i2e*^%0,t2 + he'^^^h + he-'^^y]; [ti + 
he'^yM + ^le-*^^0,^l + i2e-'^-]-[t2 + iie''^y,h + 
t2e^^y,ti + i2e''^°=,0]. For simplicity, we assume ^2/^1 — 
^2/^1 = oi, ti/ti = 12/ 12 = /3 in the following discussions 
(setting ti — 1), the extensions to more generic cases are 
straightforward. Hence, a, (3 measure the strengths of hop- 
ping frustration (due to the Se-bridged hybridization) and the 
structure distortion (due to the Fe-vacancies) respectively. The 
obtained four bands £i,2,3,4(-f'^) are illustrated via the band 
structures and the density of state (DOS) for several sets of 
(a, (3), shown in Fig. 2(a-d). The sub-band gaps Ai2,A34 
as well as the mid-band gap A23 are plotted in Fig.3(a-c). We 
find that the system is always metallic at half-filling except for 
/3 < 0.18, where the mid-band gap A23 opens at the Fermi en- 
ergy for a > 1. 

On the other hand, in the large U limit, the n.n. and 
n.n.n. hoppings lead to the AFM exchanges (Ji = Atl/U, 
J2 = Atl/U) and (Ji = Atl/U, J2 = ^il/U) for the intra- 
and inter-block spins, respectively. The resulting spin Hamil- 
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(a)a = 0.0, p = 1.0 
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(b)a=1.0, p = 1.0 
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FIG. 2: (color online). Non-interacting band structures and DOS 

p{E). 



(a) 




FIG. 3: (color online). Various non-interacting band gaps : (a) sub- 
band gap A12; (b) mid-band gap A23; (c) sub-band gap A34. 



tonian reads 



where 



(2) 



Hp = { J2 E J2}Sxip)-Sy{p) 

Hp.., = { E -^1+ E J2}Sx{p)-Sy{q) 



are the intra- and inter-block Hamiltonians, respectively, with 
p,q = (m, n) being the n.n. blocks. The spin Hamiltonian is 
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well defined on the new bipartite lattice with AFM exchanges 
connecting the even and odd sublattices. It can be further 
shown that the lowest energy state has positive definite co- 
efficients in the rotated Ising basis (obeying the Marshall's 
sign rule I ItI ITsIl ). hence the ground state is a singlet of the 
total spin Y,p Sbiock{p), where Suock{p) = ^>'(p) the 
total spin in the p-th block. This implies that for each fixed 
A-site, {S\(j)),\/p\ forms a spin lattice of period -\/5 x a/S, 
resulting in the Neel AFM orders in the classical limit [ex- 
cept for Sbiock{p) — where the ferromagnetic order is pos- 
sible Iil8i] 1. while the relative directions among different spins 
{S\{p),\/X} in a given block are dependent on the intra- and 
inter-block couplings. Notice that the block spin Suockip) 
should not be understood as a single entity, since the inter- 
action cannot be rewritten in terms of Sbiock{p)- Rather, its 
formation is protected by the periodicity of the AFM bipar- 
tite lattice and a strong intra-block coupling is not a prereq- 
uisite. The expectation value of the block spin {Sli^^j^{p)) 
varies in between and 2 and could reach the maximal mag- 
nitude when 2 J2 + Ji > 2Ji + J2. 

We now consider the weak and intermediate U and solve 
the model by the MF method as applied to the conventional ti- 
t2 Hubbard model [l^iil. We introduce the magnetic order 
parameters Mx by writing 

{2Sl{p)) = Mp)) - Mp)) = Mx cos(Q • r{p)) (3) 

for each A-sites, with r{p) = (m, n) being the position of 
the p-th block. We also consider the wavevector Q = (tt, vr), 
i.e., the Neel AFM ordering, for there is no direct coupling 
between the n.n.n. blocks. In our model, the coUinear AFM 
ordering Q = (tt, 0) (or Q = (0, vr)) may be favored when 
i^biockip)) ~ '^^'^ this case is simply identified with the 
collinear AFM in the ^2-^2 Hubbard ( or J1-J2 Heisenberg ) 
model (without vacancies) favored for strong frustration t2 (or 
J2)|l2il231. 

The MF approximation for the on-site repulsion U -term is 
standard, and diagonalization of the MF Hamiltonian is im- 
plemented over the first magnetic BZ zone of A'-space (the 
reduced BZ for one of the sublattice). Eight bands of quasi- 
particles Ei^2,---,8{K) are obtained, and the magnetic mo- 
ments Mx as functions of U are self-consistently determined 
by tuning the chemical potential ^ at the half-filling. We 
first focus on the moderate frustration regime (a ~ 0.8) with 
weak, moderate, and strong inter-block hopping amplitudes 
exempUfied with (a)/3 = 0.2, (b) /3 = 1, and (c) /3 = 1.5, 
respectively. For each case, the magnetizations Mx as well as 
the magnitude of the mid-band gap A at the Fermi energy as 
functions of U are shown in Fig. 4(a-c). We find that the four 
Mx=A,B,c.D are equal to each other in the calculated cases 
whereas they could be different in more generic cases. To 
draw the phase diagram, we perform the same calculations by 
varying the frustration a ^ 0—1.5 for each fixed /3. The AFM 
instability Um and the MI transition point Umi are determined 
by the onsets of magnetization and gap within the accuracy 
~ 10^^. For comparison, another point L/vh is marked at the 



kink position in the magnetization curves in Fig. 4. This point 
is identified as the van Hove singularity at the Fermi surface as 
evidenced by the corresponding band structures(not shown). 
Then, the a ~ U phase diagrams are plotted in Fig. 5 for each 
cases. Notice that in Fig. 5(c) the AFM metal region above the 
dashed line (a < 0.8 ) may actually merge into the PM phase 
because the magnetization is vanishingly small within the cal- 
culation accuracy. Finally, three main aspects of these results 
are briefly discussed as follows. 

(I) For moderate inter-block hopping {(3 ^ 1), the ground 
state is a PM metal for sufficiently small U and an AFM in- 
sulator for relatively large U. An extended intermediate AFM 
metallic phase emerges, separating the PM metal and AFM in- 
sulator phases in a wide range of frustration a ~ 0.5—1.5. For 
instance, when a — 0.8, we find Umi ^ 4.4, and Um ~ 3.0, 
respectively, while the van Hove singularity, which is located 
in between Um and Umi, could be clearly identified. It is 
worth noting that in the corresponding fi-<2 Hubbard model 
without vacancies the Neel AFM metal phase is extremely 
narrow and is replaced by the collinear AFM metallic phase 
for a > 0.7|22]. The corresponding Umi ^ 5.0, Um ~ 3.6. 
So we readily find that the Coulomb U needed for the Mott- 
transition is lowered from Umi ~ 5.0 to Umi ^ 4.4. A sim- 
ilar conclusion for Umi is obtained by the spin-rotor MF the- 
ory for the two-orbital model appropriate for the iron-deficient 
ternary iron chalcogenides[9j. 

(II) In the strong inter-block hopping regime /3 ^ 1.5, the 
distinction between the PM and AFM metal phases is not sig- 
nificant, in particular for a < 0.8, where the magnetization 
is vanishingly small. Thus in this regime, the AFM ordering 
immediately leads to the MI transition, and the van Hove sin- 
gularity f/vH is also very close to Umi- Another feature is that 
the Umi shows a minimum in the moderate frustration regime 
around a ^ 0.8, while the gap dependence on U is almost 
linear after the minimum. This feature resembles to that in the 
moderate inter-block hopping regime, see in Fig.4(b-c). 

(III) When the non-interacting system is close to the band 
insulator phase, as in the case of weak inter-block hopping 
P ^ 0.2 and for a < 0.8, we find that the band-gap opens 
at the Fermi energy immediately with very small Umi, after 
the development of the Neel AFM order The magnetization 
Mx rapidly saturates to 1/2, and the gap magnitude A, while 
being very small, keeps almost unchanged for further increase 
of U below a threshold value U ~ 4.5, manifesting itself as 
a conventional AFM semi-conductor. However, a drastic in- 
crease in Mx is observed above this point, where it rapidly 
saturates further to the maximal value (defined as unit), see in 
Fig. 4(a). Meanwhile, the gap increases almost linearly with 
the same tendency seen in (I) and (II). Moreover, the PM 
phase shrinks and Umi increases for a > 0.8, so an interme- 
diate AFM metal phase develops substantially in the strong 
frustration regime. 

To summarize, we study the single-orbital Hubbard model 
on the square lattice with the tetragonal vacancy super- 
structure, the most perfect vacancy ordering pattern ob- 
served in the recently discovered iron deficient compounds 
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FIG. 4: (color online). ?7-dependences of magnetic moments Mx and the mid-band gap A. 




FIG. 5: (color online). Phase diagrams. The yellow, red, and green regions indicate the PM metal, AFM metal, and AFM insulator phases, 
respectively. Notice that in (c) the magnetization in the AFM metal region (pink region(light gray)) above the dashed line for a < 0.8 is 
vanishingly small so that this AFM metal region may actually merge into the PM region (yellow region). 



(A,Tl)yFe2-a:Se2. With increasing Coulomb interaction U 
the model shows a block-spin type Neel AFM order as a con- 
sequence of the vacancy symmetry. The MI transition takes 
place at a relatively small Uui compared to the conventional 
Hubbard model without vacancies. But in a wide range of pa- 
rameters an extended intermediate region of the metal phase 
with block-spin Neel AFM exists, separating the PM metal 
and AFM insulator phases. This feature is also in contrast to 
the models without vacancies. In some cases, an almost linear 
gap dependence on U accompanied by a broad crossover of 
the magnetization is observed slightly above Umi and this be- 
havior is indicative of the electron correlation effect within our 
MF approximation. An interesting exception is the case for 
the weak inter-block hopping regime, where drastic increases 
in both magnetization and gap magnitude are observed in the 
insulator phase. Previous LDA calculations suggest that inter- 
block n.n.n. AFM exchange J2 dominates in the ground state 
of (A.TI)n ><Fei fiSe^ lillil . indicating that the system should be 
in the strong inter-block hopping regime. While the multi- 
orbital characteristic and the Hund's rule coupling should be 
also taken into account for realistic materials, our results shed 
new Ught in understanding the electron correlation and the 
novel magnetism in the (A,TI)j,Fe2-2;Se2 compounds. 
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